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Abstract. Using a reweighting technique combined with a low- mode truncation of the 
fermionic determinant, we estimate the quark-mass dependence of the QCD topological 
susceptibility with overlap fermions. In contrast to previous lattice simulations which all 
used non-chiral fermions, our results appear to be consistent with the simple continuum 
model of Diirr. This indicates that at current lattice spacings the use of the index theorem 
might not be justified and the fermionic definition of the charge might be needed. 

Smooth SU{3) gauge field configurations can be characterised with an 
integer charge Q, a topological invariant. Topology plays a crucial role 
in shaping the low-energy behavior of QCD. Instantons are believed to 
be largely responsible for the spontaneous breaking of chiral symmetry. 
The axial anomaly [||] as well as the large mass of the rj' ^ are also 
intimately connected to gauge configurations of non-trivial topology. For a 
good understanding of QCD, it is thus crucial to have a consistent picture 
of its topological structure. 

Fluctuations of the topological charge can be characterised with the 
topological susceptibility, 

x^^. (1) 

where (.) denotes averaging with respect to the full QCD measure and 
V is the volume of the system. One of the most profound effects that 
light dynamical fermions are expected to have on the QCD vacuum is the 
suppression of fluctuations of the topological charge. Chiral perturbation 
theory predicts that in the presence of Nf degenerate light fermion flavours 
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the susceptibility vanishes with the quark mass as 

where and are the pion mass and decay constant, S is the chiral 
condensate and m is the quark mass Q. 

At the other extreme, when the quarks are very heavy, their influ- 
ence eventuahy becomes neghgible and the susceptibihty approaches its 
quenched value which is known to be Xq = (203 ± 5MeV)^ |5[. Recently 
it has also been discussed by S. Diirr how the susceptibility is expected to 
behave for intermediate quark masses Q . Starting with the observation that 
besides the light quarks, the unit volume^ also suppresses higher topological 
sectors and noting that these two mechanisms are independent, he derived 
a phenomenological formula for the susceptibility. In terms of the pion mass 
it reads as 

X{m^)=(^% + —] ■ (3) 

Xq) 

Both and Xq &re independently known therefore this formula has no free 
parameter. While ^ does not come from flrst principles, it is based on a 
physically appealing plausible picture. Even if it is not the last word in this 
subject, it indicates that the susceptibility might substantially deviate from 
the Leutwyler-Smilga curve already at smaller quark masses than previously 
expected and it might be strongly suppressed even at moderately heavy 
quark masses. 

Although it is only the quenched topological susceptibility that has 
an immediate phenomenological interest, being connected to the r]' mass 
through the Witten-Veneziano formula [§, ^, it would also be of consider- 
able interest to check the expected chiral behaviour of the full (unquenched) 
susceptibility. This is because in the derivation of the Witten-Veneziano 
formula it is always tacitly assumed that the full susceptibility goes to zero 
in the chiral limit (see e.g. for a recent discussion on the lattice). 

Recently several attempts have been made to check these predictions 
against numerical simulations on the lattice but the situation is still rather 
controversial. Some lattice simulations show only very slight or no suppres- 
sion of topological fluctuations |8|, 10 1. The latest UKQCD results [11| 



obtained with improved Wilson fermions agree well with the Leutwyler- 
Smilga prediction of Eq. (Q), even beyond its expected range of validity 
and the UKQCD susceptibility is still signiflcantly higher than what Diirr's 
interpolation formula (^) anticipates. 



^ Recall that the susceptibility is defined as the fiuctuation of Q per unit volume. 
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To understand why lattice simulations do not show the expected fermionic 
suppression of the topological susceptibility, let us briefly look at the physi- 
cal mechanism that leads to this suppression. In the background of a charge 
Q smooth gauge configuration the Dirac operator D has (at least) \Q\ zero 
eigenvalues for each flavour, due to the Atiyah-Singer index theorem. The 
quark effective action — obtained after integrating out the fermions — for 
Nf flavours is proportional to det^^(L' + m). In the light quark limit the 
determinant is oc m}^^^f which leads to the suppression of higher Q sectors 
and thus the topological susceptibility. All this is true for smooth gauge 
fleld configurations in the continuum. The situation in lattice simulations 
however is rather different since lattice Dirac operators in general do not 
have exact zero eigenvalues and it is not obvious how the index theorem is 
realised on the lattice. 

For a more thorough understanding of what happens on the lattice, 
we first note that both the Witten-Veneziano formula and the Leutwyler- 
Smilga relation follow from the flavour singlet axial Ward Identity (WI). A 
general (Wilson type) lattice fermion action can be written as 

S{e) = Tl)^ + e^'^'-'^iW + rricM) + (m - mc)Mj V, (4) 

where is the naive lattice Dirac operator, W is the Wilson term, M 
the mass term and we already split off the critical mass rUc to account 
for the non-trivial mixing between the Wilson and the mass term, and 
the ^-dependence is also included |12]. Using that anticommutes, while 



the rest of the Dirac operator commutes with 75, it is easy to show that 
the ^-dependence can be transferred to the mass term by a chiral rotation 
V'(x) e-''^''^/^i;{x). This also means that in the chiral limit, m mr, 
the action does not depend on 6. 

The flavour singlet axial WI can be derived by making a local change 
of variables ip{x) — > e~'^"'^'^^^^ilj{x) in the fermionic path integral and using 
that the measure is invariant with respect to this. The variation of the 
different terms in Eq. (Q) give rise to the following terms in the WI, 

a^j^(x) - 2Nf q{x) + 2(m - m,)Nfi^j5tl;{x) = (5) 

By integrating this over all space-time, the divergence of the axial current 
does not give any contribution and we arrive at 

Q = — (m — mc)Tr ^5D{m)^^ , (6) 

where Q is the integral of q{x) over space-time. The term Q that originates 
from the Wilson term in the action, is usually identified as the topological 
charge [O]. We would like to emphasise, however that this identification 
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is based on perturbative arguments that are valid only in the naive con- 
tinuum limit and might also hold in the proper continuum limit but there 
is no reason to believe that any charge obtained directly from gauge field 
observables (e.g. by cooling) is a reasonably good approximation to Q in 
practical situations. 




Field theoretic 

Figure 1. The value of the field theoretic vs. the fermionic charge on a set of (3 - 
quenched gauge configurations. 



5.85 



In fact, there is good reason to believe that at current values of the 
lattice spacing, these quantities are quite different. In Fig. |l] we show a 
particular field theoretic charge |5| versus Q on a set of /3 = 5.85 quenched 
configurations!^. Even though, there is a strong correlation between Q and 
the algebraic charge, Q, it is not unlikely that the use of Q instead of the 
"fermionic" definition Q can introduce substantial errors. 

A possible way to avoid this would be to use the fermionic definition, 
with the Dirac operator that enters the dynamical simulation. There is 
even an efficient algorithm to calculate Q [0] but difficulties emerging in 
the chiral limit have so far made this approach impracticable. 

Recently, through the overlap formulation, it has become possible to 
put fermions on the lattice while preserving an exact chiral symmetry [15| 
and avoiding many of the problems other formulations have in the chiral 
limit. In particular, with the overlap, the above derivation of the WI and 
the fermionic definition of the charge results in 



n_ 



n+, 



(7) 



where n± are the number of positive/negative chirality zero modes of D. 
This makes it possible to define the "topological charge" in a consistent 



^ The Q values have been slightly rescaled by an overall factor to give the same 
topological susceptibility as the field theoretic charge. 
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way avoiding any reference to an algebraic charge that has no a priori 
connection to the Dirac operator. Another advantage of this formulation 
is that with chiral fermions the Leutwyler-Smilga relation has been shown 
to hold exactly in the chiral limit even at finite lattice spacing This 
makes chiral fermions extremely suitable for exploring how the topologi- 
cal susceptibility interpolates between the chiral Leutwyler-Smilga and the 
heavy quark regime. 

Unfortunately lattice simulations with dynamical overlap fermions have 
so far been prohibitively expensive. In the present paper we show the results 
of a qualitative numerical computation of the susceptibility with overlap 
fermions. As we shall see, the quark-mass dependence of the susceptibil- 
ity can be estimated with a suitable reweighting of a quenched ensemble. 
The procedure we propose is to generate an ensemble {Ui}i of quenched 
gauge field configurations with only the pure gauge measure. In principle 
the topological susceptibility can then be computed by reweighting each 
configuration with the corresponding fermion determinant, 

Xi^) = ^ E 2' det^/ [Dm + m], (8) 

i 

where D{Ui) is the Dirac operator and Qi is the charge in the background 
of the gauge configuration C/j. Z is the partition function, i.e. the sum 
appearing in Eq. (^) without the factor. 

In principle this is a correct procedure, in practice however it is useless 
since control over the statistical errors is lost exponentially with increasing 
volume. The crucial observation is that it is only the small eigenvalues of 
D that are correlated with the topological charge and without any loss, the 
full determinant can be replaced with a truncation thereof, 

N 

det(i? + m) ^ [](Afc + m), (9) 

fc=i 

where {Ai, A2, • • • Aat} are the N smallest magnitude eigenvalues of D. 

In Fig. |2| we show the (reweighted) susceptibility as a function of the 
number of eigenvalues included in the truncation. For heavier quarks, more 
eigenvalues should be included in the determinant but overall the trunca- 
tion is clearly a good approximation. It drastically reduces the variation of 
the determinant within each topological sector and turns out to make the 
reweighting practicable in a physically useful range of volumes. The upshot 
is that once the lowest N eigenvalues of D are known on an ensemble of 
quenched gauge configurations, the full quark-mass and Nf dependence 
of the susceptibility can be easily obtained. The only caveat is that for 
smaller quark mass and larger Nj the determinant fluctuates more and 
consequently the errors increase. 
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Figure 2. The reweighted topological susceptibility in a = 2.29 fm* box, as a function 
of the number of eigenvalues included in the truncated determinant. The different curves 
correspond to different quark (and pion) masses. 



We have seen that the truncation of the fermion determinant is a good 
approximation for the susceptibihty but what does it physicaUy mean to 
truncate the determinant? We would hke to argue that it is essentially just 
another possible action that differs from the untruncated fermion action 
only by local gauge terms. Indeed, it has been shown that the contribu- 
tion of the upper part of the spectrum to the fermion action can be well 
approximated with a linear combination of small Wilson loops [17|, i.e. a 
local gauge action. 

Besides the susceptibility, there are two other quantities of interest in 
our simulation, the lattice spacing (as set by the Sommer scale ro) and the 
dependence of the pion mass on the bare quark mass. It is expected that 
including the truncated determinant will not substantially change tq com- 
pared to its quenched value, in conventional dynamical simulations most 
of the change in tq originates from the bulk of the upper part of the Dirac 
spectrum. Therefore, in this qualitative calculation we completely ignore 
this effect and use the quenched scale. The change in mT^{mq) compared 
to the quenched one can be estimated using the GMOR and the Banks- 
Casher relation together with the available low-end of the Dirac spectrum. 
The setting of the scale and the pion mass are the two weak points of this 
calculation which might contain uncontrolled effects. This is the reason why 
this computation should be considered only as a qualitative first estimate 
and hopefully in the future it will be checked by a full-blown dynamical 
overlap calculation. We believe, however that even with these caveats, our 
calculation is still "competitive" with the currently available dynamical 
calculations that might have other limitations, as we discussed above. 

Simulations at several (3 values and lattice sizes are currently underway. 
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As an illustration we show results on an ensemble of 2500 /3 = 5.85 (a = 
0.123 fm) lattices with a volume of ^ = 2.29 fm*^. Our overlap action was 
constructed by inserting the Cgw = 1-0 clover operator with 10 times APE 
smeared gauge links [p!8| into the overlap formula, 



[l-Dr 



[l-D,){l-D,)^ 



(10) 



The inverse square root was approximated with Chebyshev polynomials 
after treating the lowest eight eigenmodes of (1 — Dc)(l — L^c)^ exactly [19|. 
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Figure 3. The topological susceptibility vs. pion mass in a box volume of = 2.29 
fm*. 



In Fig. ^ we present the main result of this paper, the susceptibility 
versus the pion mass with Nf = 2 flavours of overlap fermions. Both the 
susceptibility and the pion mass are translated into dimensionless units with 
the help of the Sommer scale rg = 0.49 fm pO|] . The shaded region indicates 
the quenched susceptibility r^Xq = 0.0703(19). obtained on the same ensem- 
ble without any reweighting. The quoted error is only statistical, we did not 
include any additional uncertainty coming from the scale setting. In more 
conventional physical units this translates into Xq = (207.0(1.4)MeV)^. 

For comparison we also plotted the Leutwyler-Smilga formula and Diirr's 
formula of Eq. (^ using our measured quenched susceptibility and the 
known value /^r = 93 MeV. Since both Eq. (|3|) and (|2|) are valid in the 
infinite volume limit and our volume is not particularly big, we also plotted 
the curves corresponding to our volume. 

In conclusion, we argued that current computations of the topological 
susceptibility with dynamical quarks might suffer from a systematic error 
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essentially due to the lack of the index theorem at currently used lattice 
spacings. We presented a pilot study of the unquenched topological sus- 
ceptibility with Nf = 2 flavours of overlap fermions using reweighting. 
Our computation involved some uncontrolled approximations for setting 
the scale and the pion mass which eventually have to be checked in a 
full-fledged dynamical simulation. Nevertheless, our results appear to be 
consistent with Diirr's formula and indicate that the susceptibility might be 
more strongly suppressed than other lattice simulations suggest. It would 
also be interesting to compute the susceptibility with non-chiral Wilson 
fermions using the corresponding fermionic definition of the charge. 
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